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On the Problem of Two Fixed Centres and Certain of 

Its Generalizations. 

By Adam Miller Hiltebeitel. 



The classic problem of two fixed centres may be stated thus: To determine 
the motion of a body attracted by two fixed Newtonian centres of force. It was 
first solved by Euler and since Euler generalizations of the problem have been 
made. In this paper are given the resolution of the most general two-centre 
problem of which the variables are separable and a discussion of the trajectories 
of a less general form of the problem. 

Problem : To study the motion of a material body M acted upon by five 
centres of force, K lf K % , K 3i K it K By which are the two real foci, the centre, and 
the two imaginary foci, respectively, of a system of confocal conies, by forces 

R l ——mr 1 — -J, B z = — m,r s —-^, B 3 = — rn 3 r 3 , 



2 



7?— m 'r m i + im z R — m i r m J— im 6 
Mi rn r 4 -g , xi h — — m r 6 -, , 

*4 '6 

respectively, and by an additional force always parallel to the line K± K h and 
varying inversely as the cube of the distance of the moving body M from the 
fixed plane perpendicular to K± K 6 at K 3 , the centres K 3 , K if K 5 being fixed, 
while the real foci K lt K 2 are always in the meridian plane K± MK 6 . 

Let x, y lf % be the coordinates of the moving body M referred to a fixed 
system of rectangular axes whose origin is at K 3 and whose tc-axis coincides with 
the line passing through the imaginary foci K it K 5 , as in Fig. 1. Imagine a 
plane through the cc-axis and the moving body if, and let the intersection of this 
moving plane with the fixed y 1 %-plane be YY'. This line will always lie in the 
meridian plane E i MK h and, as this plane rotates about the ai-axis, YY' will 
rotate about the origin. Every point in this line will describe a circle in the 
Vi %-plane, whose centre is K 3 , and upon it the two real foci are taken. Of the 
five centres of force, therefore, K 3 , K if K h are fixed, while the other two move 
on a circle whose centre is K 3 . 
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The motion of the meridian plane, and therefore that of the centres K lt 7T 2 , 
is due to that element of the initial velocity of M which is perpendicular to the 
initial position of the meridian plane. But this element of the initial velocity 
not only causes the rotation of the meridian plane, but also influences the nature 
of the path described by the moving body in this plane. By the aid of the 
principle of areas the motion of M in the rotating plane can be completely 
separated from the motion of rotation about the cc-axis, as Jacobi* has done 
for the simpler problem. 

The perpendicular distances of M from the a-axis and the line YY' are 
designated by y and x respectively. Regarding XX' and YY' as a system of 
rectangular axes in the meridian plane, then x, y are the rectangular coordinates 
of M defining its position in this plane. If at the beginning of motion the plane 
KiMK 6 coincides with the fixed plane XK S Y 1} and during the time t describes 
the angle $, as indicated in Fig. 1, then 

y x = y cos $, z 1 — y sin q>, $ + % = y*. (1) 




Fig. 1. 

Let 2c denote the distance K x K % and lie the distance K^K^. Let the 
rectangular coordinates of K lf K z in the Fj Z r plane be ( — a, — 6), (a, I), 
respectively, and the distances K x M, K% M, K s M, K± M, K & M be respectively 
r i) r 2> r 3> r i> r 6' We now have : 

a* -f- & 3 = c 3 , by 1 — az 1 — 0, (2) 

r\ = x^ + (y 1 + af+{z 1 + hf, r\- x* + {y x -a? + (%-6) 3 , (3) 

ri = x" + yl + zl rl=(x + icf + y\ + zl, r\= (x- icf + y\ + z\. (4) 

* Vorlesungen uber Dynamik, No. 29. Werke, Supplement, pp. 321-331. 
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The form of the forces R lt R 2) R if R 5 is accounted for by supposing that 
there are two masses at each of the corresponding centres, one acting inversely 
as the square of the distance and the other directly as the distance. At K lt K% 
the masses acting inversely as the squares of the distances, denoted by m lf m z 
respectively, may be any real quantities whatever, while the masses acting 
directly as the distances must be real and equal to each other. These latter 
masses are denoted by m. But at K it K 6 the masses acting inversely as the 
squares of the distances will be denoted by the conjugate complex quantities, 
<m± + im 5 , tn^ — im 6t i being the imaginary unit, and the masses acting as the 
direct distances may be any real quantity, as m'. Under these hypotheses as to 
the masses, the potential function will be real and the separation of the variables 
in the Hamilton-Jacobi partial differential equation can be effected by the aid 
of elliptic coordinates. 

This problem includes, as particular cases, all the various forms of the two- 
centre problem. The form of the potential function given by Darboux* may be 
obtained by setting m = m! = 0, or by compounding the two forces — m r x , 
— mr s into the single force — 2m r s and suppressing the forces — m'r i} — m'r 5 . 

These four forces together with the force — m s r s give rise to the following 
terms in the potential function : 

— | <m r\ — \ m r% — \ tn 3 r% — \ m' i\ — | in' r\ ; 
and it can easily be shown that this expression is equivalent to 

— (rn -f- \ m s + in') r% — c 3 {rn — in'). 

Hence, even on retaining these forces, the potential function of this problem 
differs from that of Darboux's problem only by a constant term. Velde'sf 
problem is obtained from the present problem by omitting the forces R i} R 5 
as well as that parallel to the sc-axis, and by regarding the motion as in a fixed 
plane through the centres K lt K 2 , which are also supposed fixed. By simply 
omitting the forces R if R 5 , we have the most general form of the two-centre 

* Archives Merlandaises des Sciences, Ser. 3, Vol. VI, pp. 371-376. 
+ Programm der ersten hoheren Burgerschule zu Berlin, 1889. 
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problem as formulated by Liouville.* Retaining only the forces M lt B 2 and 
regarding the corresponding centres as fixed, we have the problem as extended 
by Lagrange ;f and if, in addition, we suppress the terms — mr x , — vnr 2 , in 
R x , B 2 respectively, we have the original two-centre problem of Euler. J 
Differentiating equations (2), (3), (4), we have 



a da + b db = 0, b dy x + y x db — a dz x — z x da = 0, (5) 

;}w 



r x dr x = x dx + (y x -f- a) dy x + (y x + «) da + (z x + b) dz x + (z x + b) db, 
r 2 dr 2 = x dx + (y x — a) dy x — (y x — a) da + (z x — b) dz x — (z x — b) db 



r 3 dr s = xdx + y x dy x + z, dz x , 

r± dr 4 = (x + ic) a*.r + ^ dy x -f % c?z x , 

r 6 aV 6 = (a; — ic) dx + 3^ a*^ + z x dz x . 

But 

(y x + a)da -\- (z x -\-b)db — -±- (ada + bdb) + aoa + bdb = 0, 

— (2/1 — a) aa — (zj — b) db = — ^- (ada + &a7>) + ada + 6df6 = ; 

and therefore the equations (6) become 

r x dr x = xdx + (y x + a) dy x + (z a + &) dz x , 
r 2 dr 2 = xdx -f (y x — a) dy x + (z x — 5) dz x . 



(7) 



(8) 



It will now be necessary to express the distances in terms of the variables 
x, y. This can be done without difficulty by aid of the relations (1) and (2), 
and we have 

r\ = x 2 + (y + c)\ ri = x 2 + (y — c) 2 , r\ = x 2 + y 2 , "| 

r\ — (x + ic) 2 + y 2 , rf = (x — ic) 2 + y 2 . J ( ; 

* Joum. de Math., Vol. XI, pp. 345-378; Vol. XII, pp. 410-444; Vol. XIII, pp. 84-37. Oonnaismnce de 
Temps pour 1849, pp. 255-256. 

I Miscellanea Taurinensia, Vol. IV, pp. 118-243 : (Euvres, Vol. II, pp. 67-121 ; Micanique Analylique, 1st Ed., 
pp. 262-286; 2nd Ed., Vol. II, pp. 108-121; (Euvres, Vol. XII, pp. 101-114. 

X Nov. Comm. Acad. Imp. Petropolitanae, Vol. X, pp. 207-242; Vol. XI, pp. 152-184. MSmoires de V Acad, 
de Berlin, Vol. XI, pp. 228-249. 
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Taking the mass M as the unit mass, the differential equations of motion 
may now be written as follows : 

d z x is , • \ it • \ wJia; m 2 x 

T-ir = — mx — mx — m % x — m (as 4- tc) — m (x — ic) ^ 4- 

dv v ' v ' r\ rl 

(ot 4 + im 6 ) (x + id) (ot 4 — irn 6 ) (x — ic) {3 



X s ' 



V± = — m(y x + a) — m{ yi — a) ~ ni 3 y 1 — m'y x — m'y x 



dt 



m x (y x + a) _ m z {y x — a) _ {m i + im^y x _ (wt 4 — im s )2/i 



rl r\ r\ 



-j-±- = — m (z x + b) — m(z 1 — b) — m 3 z x — m! z 1 — rn' z x 

_ m 1 (z 1 + b) _ tn z (z x — b) _ (rjij + im 5 ) z± _ (w 4 — im 5 ) z x 

rl rg rl "' r 6 8 

Multiplying equations (10) respectively by dx, dy x , dz x , adding and inte- 
grating, we obtain the vis viva integral 



(10) 



r x r z r 4 r 6 2 a? 



(11) 



Now multiplying the second equation of (10) by z X) and the third by y x , 
and then subtracting, we obtain 



_ "hj%i — agfi) i ™»{byi — az x ) . 
and therefore, since %! — az a = 0, we have 



(12) 



y d 2 z x _ d^y x _ 

yi ~W Zl ~dF~°- ( 13) 

Integrating this equation, we have 

dz x dy x ,, jN 

where a is the constant of integration. This integral equation shows that the 
areas described by y, as the meridian plane rotates about the a>axis, are pro- 
portional to the times. 
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By the aid of (1) we find 

(t) ! +(l') ! =(l)'+^@)'. (»> 

Substituting these values in (11) and (14), we have 



, vn 1 . m 2 m i + itn 5 , 'm i — im 5 , 1 §__ x a 2 , h 
»i »» n ~ r r 6 ' r2 'x 2 2 'y 2 ~ r ' 



(17) 



V*f = a, (18) 



in which the variables y lt z lf a, b no longer appear. 

Differentiating the equation y = s/y 2 + %\ twice with regard to t, we get 

*,_ *&+*% > , w^0^(S)|-(»l +'■§)' (,., 

dt 2 ~ y y s 

The second term of the right-hand member of this equation reduces to 

1 / dz 1 %A 2 _ a 2 

y-*V l ~di~ Zl dt) -jf' 

and the first term becomes 

— m (y + c) — »n (y — c) — w? 3 ?/ — m'y — m! y 

m i {y + q) _ ^z (y — q) («»4 + Jm s ) y Qw 4 — t«i 8 ) y 
»*i 3 »? n 3 rg 

Substituting these values in (19), we find that the original equations of motion 
may be replaced by the following system : 

d 2 x 



dt 2 



= — rnx — mx — m s x — m' (x + *c) — #?/ (x — id) 

_m 1 x rn 2 x_ (m. 4 + **%) (g + jjg) (m 4 — i?n 6 ) (a; — io) /? 

/f ri rf ~rj X s ' 

-to 

-^V- = — m(y + c) — m(y — c) — m 3 y — m'y — m'y 

m^y + c) m 2 (y — c) (m 4 + im 6 ) y (mj—im^y a 2 
r? ~~r$ r\ r 6 8 y 3 ' 

dp a 

~dt y~ 2 ' 



(20) 
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The first two of these equations involve only the variables x and y and 
define the position of the moving body in the meridian plane. The third 
equation gives the motion of rotation about the a;- axis. 

If we set 

U = — \ mrf — \mrl — \ r m 7i rl — \ m' rf — J m! rl 

rn 1 Wg (w 4 + im 6 ) , (w 4 — im h ) , x §_ _ .. i « 2 

n n n "*" r 6 " t " 2 x 2 2 y 2 ' 

the equations (20) may be written in the form 

<Px _dU d*y_dU d$__dU 
dt* " dx ' dt z dy ' dt da' 

On introducing elliptic coordinates X, (i defined by the equations 

^ = i(r 1 + r s ), (1 = 1^! — r,), (21) 

the functions T and ?7 become 

-=H^,-){^(§y-^(*y}, 

U = ^-L_ a |(- (m + | m, + m>)tf + c*{\ rn 3 + 2 m')^« + K + «,)*. 
+ 2 W , Vl^=^ + | jj^s- l°^- 2 ) + ((m + | ™ 3 + m')i« 4 
-^(im 3 + 2mO^-K- ma )^-2m 6 V^^-i (( ^+i^ 2 )}, 
by virtue of the relations 



»-! = *, + /«, r 2 = a — p, r 8 = V ^ + ii 2 



r 4 = V/l 2 — c 2 +V,k 2 — c 2 , r 8 = V 7? — c 2 — */~j? — <?, 
c z y 2 = ^(i z , c 2 o; 2 = (a 2 — c 2 )(c 2 — ^ a ). 

Both T 7 and U have the form required for the separation of the variables 
in the Hamilton-Jacobi equation. U, however, is an irrational function of the 
coordinates % and (i, but it is real for all real values of % > c and for all values 
of (i numerically less than c. It is a real function, therefore, for all values of 
% and fi consistent with their definition. 
45 
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Setting H=T — U, the differential equations defining the motion in the 
meridian plane assume the canonical form 

dX _dH dpi_^ __dH 
dt dp^ dt dX ' 

dfx dH dp 2 327 

dt ~~ dp z ' dt ~ d(i ' 

where the conjugate variables^ and p 2 are given by the relations 

_ dT _dT 

Pl ~dX" Pz ~d^' 

The Hamilton-Jacobi equation, 

i f a 2 — c' fdW\ a u 2 — c 2 /dW\ z \_ TT ,j i 

becomes, on inserting the value of JJ in full and arranging, 

x (# _ ») (|^) 2 — i— (m + I n> 3 + m')X + (± m 3 c* + 2 m'c 2 + h) X 2 

1 _/3c a _ _ia^l 
2 X 2 J 



+ (wi! + Wl a ) X + 2 OT 4 V ^ 2 — c 3 + 

'3TTV 



2 A 3 — c 2 2 



-(m 1 - Wl2 )^-2 W6 V^^-|^^ + J^£ 2 

In this equation the variables are separated. Hence by setting each member 
equal to the same constant and putting for short 



A=/(^) + (»i 1 + m 2 )a 5 + 2m 4 aV^-c 2 -2m i c 2 XV^-c 2 

— (m 1 -{-m 2 )c 2 X z , 



M.=f{n) + {nii — m2)ti 5 —2m 6 nWc z —ii z + 2m 6 c 2 iiWc z — [i !> 

— (wjj — wi 2 )c 3 ^ 3 , 

/(«) = — (m + i »i 3 + m') s 8 + (wc 3 + m 3 c 2 + 3 m'c 3 + £)s 6 

— (i m 3 c 4 + 2m'c i + hc 2 —k)s i — (k + | a 3 — J/8) c 3 s 3 + a a 3 c 4 , 

we obtain the two equations 

dW V~2 ,— dW _. V~2 



a a, a(x 2 — c 2 ) 

to determine the function TF. 



dfi (i ((i z — c 2 ) 



VM 



(22) 
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Multiplying equations (22) by d\ dfi respectively and integrating, we find 

W - W(X) + W{ji). 
Two integrals of the equations of motion then are 

' + 7l _ viV VT + ;t?./ 71' (23) 

^^VA + VlivM (24) 

the first of which defines the time corresponding to any position of the moving 
body in its path and the second the path described in the meridian plane. 

The intermediary integrals (22) define the velocities, and they may be 
written as follows: 

(* - ^ § = ^ vx (Jll _ ^ dji = _ ^2 vir> (25) 

From (23) and (24), we find by differentiation 

V2dt = ^ + ^M, (26) 

VA T VM y 

VA + V M ~ °' (27) 

and with the aid of (26) and (27) we derive from the third of the equations (20) 
the relation 

V V2UVA nVMS' k ' 

and after integration 

^ +ai =_^iJ/*_^ =+ f-^j, (29) 

T V2 1V WA •/ ^VMj' 

which defines $, the angle of rotation of the meridian plane. With this integral 
the solution of the problem by quadratures is complete. The three final integrals 
are (23), (24), (29), and the six arbitrary constants of integration are h,h,h u k lt 
a, %, their values depending upon the initial conditions of the motion. 



d( 
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By their definition (21), the coordinates \ (i must always satisfy the in- 
equalities 

oo ^>%>c, — C>[l>C, 

and in order that there may be real motion, they must also satisfy the conditions 

A>0, M>0, 

not only at the beginning of motion but also throughout the motion. Further, 
X can never take the value + c nor approach + c indefinitely, for in this case 
the potential function U becomes infinite. For the same reason, (i can never 
take the values ±c nor zero, nor approach these values indefinitely. Therefore 
/I must have a minimum value, say a lt greater than + c, and ^ must have a 

Y 




Fm. 2. 



maximum and a minimum value in each of the intervals (0, -f- c), ( — c, 0). Let 
the maximum and minimum values in the interval (0, + c ) be respectively 
b x and b 2 , and in the interval ( — Q» 0), fix and /? 2 respectively. Again, if we 
suppose all the masses to be essentially positive, then for increasing values of X, 
the function A finally becomes negative and remains negative. Hence /I can not 
increase indefinitely and the maximum value of X for which A is positive or zero 
we shall denote by a 2 . Hence X and ft must satisfy the following conditions : 

a % >X>a lt hi > /x > b z , /?! > /m > /? 2 . 

The path described by the moving body if in the meridian plane, therefore, 
must always remain within one of the four trapezoidal regions of this plane 
bounded by ellipses X = a 1} X = a z and by the hyperbolas /x = b 1} [i = b %} or by 
the hyperbolas \i = (3 lt (* = (3 S , as indicated in Fig. 2. 
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If, at the beginning of motion, M is in one of these four regions, eay I, 
it must always remain within this region. As the meridian plane describes 
the angle <p, the region I will generate a tube-like region in space and within 
this space the motion of M must take place. Whatever the nature of the 
trajectory within this space may be, it can never extend beyond it. 

If 6j = ^ 2 , 6 2 = ^j numerically, then as the meridian plane rotates about 
its axis, the spaces generated by the regions I and II will be the same as those 
generated by IV and III respectively, but if b x =f= ($ % , b z zp ^ numerically, 
this will not be the case. It may happen, however, that the spaces generated 
by I and IV, or II and III, overlap ; then, supposing M to be in I at the 
beginning of motion, the trajectory is limited by the combined spaces generated 
by I and IV after a semi-revolution of the meridian plane. 

Let a and b be multiple roots of the equations A = and M = respect- 
ively. Then, if at the beginning and throughout the motion we have h = a , 
the motion will be upon an arc of the ellipse 2, = a ; and, as the meridian plane 
rotates about the sc-axis, this ellipse will generate an ellipsoid, upon a certain 
portion of which the body M must move. Likewise, if we have (i = b at the 
beginning and throughout the motion, the body M will move upon an arc of the 
hyperbola fi = b in the meridian plane, and in space, upon a certain portion of 
the hyperboloid of revolution generated by the hyperbola branch (i= b . If we 
have at the same time 2, = a and (i = b , no motion will be possible in the 
meridian plane and the orbit in space will be a circle whose axis is the a;-axis. 

A further discussion of the conditions for the kinds of motion just mentioned 
will be given later. 

When the motion is supposed to take place in a fixed plane passing through 
the five centres of force, that is, when a — 0, we have 

U -^—M~ (w» + |w 3 +m')^+(|w.3+2 m ')c 3 ^ + (m 1 + m 2 )^ 
+ ( OTl _ OTg )^_2OT s Vc a — fi 2 — $. f C )}, 



A = f{X) + {m x + m 2 )% s +2m i Vtf—d i ) 



M=f(ji) + (nti — m. z ) ,u 3 — 2 m 6 V c 8 — yp , 
/(«) = — (w + !m 3 + m')s 6 +(mc 2 + m 3 c 2 + 2ra'c 2 + A)s 4 + &s 2 +i/?c 2 — kc 2 
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When (i becomes zero, U becomes infinite. When either Xor^ becomes 
numerically equal to c, U remains finite, but again becomes infinite when 
X = (i = c ; that is, when there is a collision of M with the real focus, K x , 
or K z . Hence the trajectory may cross the y-axis but not the a-axis, and it can 
not pass through the real foci K x , K z . 

For sufficiently large values of X, the function A becomes and remains 
negative as X continues to increase. If we represent the greatest value of X 
for which A changes from positive to negative values by I, then X must satisfy 
the inequality 

+ c < X < I 
in all cases of real motion. 




Fig. 3. 

Since U becomes infinite when (i becomes zero, fi can not approach zero 
indefinitely. Representing the lower limit of the positive values and the upper 
limit of the negative values which ^ may take by 7 1 and l z respectively, (i must 
always satisfy the inequalities 

in all cases of real motion. 

From these limitations on X and {i, it follows that the trajectory must lie 
wholly within one or the other of the two domains of the fixed plane, bounded 
by the ellipse X = I and the hyperbola branches (i = l 1} (i = l z , as in Fig. 3. 

For 2, = c the ellipse degenerates into the straight line segment K x K z ; for 
fi == c the hyperbola degenerates into the straight line K x K z beyond K z , and for 
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(i = — c into the straight line K x K z beyond K x . Hence, if at the beginning 
and throughout the motion we have /I = c, the body M must move on the line 
K L K z , either on the segment B a K % or the segment B x K 1} not including the 
limiting positions K z , K x . On the other hand, if at the beginning and through- 
out the motion we have (i — + c, the body M must move upon the line segment 
K 2 A Z) excluding K 2 ; and similarly for (i = — c, the motion takes place on the 
straight line segment K x A lt excluding K x . 

If we have at the beginning and throughout the motion "k = a , where a 
is a multiple root of A = 0, the motion will be upon an arc of the ellipse /I = cr . 
Likewise, if we have [x = b at the beginning and throughout the motion, the 
motion will be upon an arc of the hyperbola (i = b . 

For a fixed system of values of the masses involved, the existence of multiple 
roots of the equations A = 0, M = depends upon the arbitrary constants h and h 
that is, upon the initial conditions of the motion. From this it is evident that 
the same conic which can be described by all the forces acting at the same time 
can also be described when any one of the forces B lf B % , B 3 acts alone or when 
the two forces B i} B 6 alone act. This fact Lagrange considered most remarkable 
and called attention to it in his work on the problem of two fixed centres in the 
MScanique Andlytique. Legendre* also noticed it. The above result and the 
theorems of Lagrange and Legendre are corollaries of the following general 
theorem due to Bonnet : f 

Si plusieurs masses m, m', m", . . . . , respectivement soumises a Paction des 
forces F, F', F", . . . ., et partant toutes d'un point A avec des vitesses v , v' 0} v'q, 
. . . ., de grandeur diff6rente mais de meme direction, decrivent la meme courbe 
ACB, la masse quelconque M, soumise a Paction de la r^sultante des forces F, 
F', F", . . . - , et partant du point A avec une vitesse V ayant la meme direction 
que les vitesses v , v' 0) v'q, . . . ■ , d^crira encore la courbe ACB f pourvu que les 
forces F, F', F", . . . . , independantes du temps et que la force vive initiale MY 
de la masse M soit egale a la somme 

des forces vives initiales des masses in, in', in", 

* TraiU defonctions elliptiques, Vol. I, p. 426. 

+ Jburn. de Math., Vol. IX, pp. 113-115 ; (Euvres de Lagrange, Vol. XII, Note IV, pp. 353-355. 
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We shall now consider the motion of a particular case of the preceding 
problem by the method employed by Charlier * in the discussion of the motion 
of the original problem of two fixed centres when the motion is in a plane 
through the line of centres. The problem is as follows : 

To determine the motion of a free body M in space attracted toward two 



Wi 



m, 



fixed centres K y , K 2 by forces B x = — m r x 1 , B z = — mr 2 f respectively, 

■ r-y r z 

and by a third force R s = — -^ parallel to the line K x K % . 




Fig. 4. 

In this case K 1} K % are fixed. The plane K X MK % is supposed to rotate 
about the line K x K % as axis, coinciding with the fixed plane XK Z Y x at the 
beginning of motion and describing the angle $ in the time t, as indicated in 
Fig. 4. AM, or x, is the distance of the moving body M from the fixed plane 
Y X K Z Z X . Hence the third force, JS 3 , varies inversely as the cube of the distance 
of M from the fixed plane T x K s Z x . 

* Mechanik dies Himmels, Vol. I, pp. 117-163. 
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If we take the mass M as the unit mass, the equations of motion are 

d z x i , \ / % mAx + c) m«(x — c) 8 

-y-g- = — m (x + c) — rn {x — c) v 8 ?A - s — K , 

dt z v ' v ! r? rl X s ' 

d z y x _ _ gh ffi __. ^a yi 

^ 2 - mfc my, ^ 3 r , , 

dz, tn-yZi m 2 y z 

dt z * 1 T-j 8 r| 

The further details of the solution are altogether analogous to those of the 
general problem already solved, so that we content ourselves with the results 

f i . / % i« 2 c 2 i/?c 2 \l 

dq> a 



w =fm^)^ + f¥i^) vw 
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Mtf-^)g= VI, (30) 

^ 2 -V)^f = -V"M, (31) 

d VdX jSdi* 

dt ~ SA + VM' (32) 

VA + Vltf- ' (33) 

^ = - 2 ((,_\yl + 0«W)Vm)' (84) 

♦ + - = --</c?^^ + /o?=^> (37) 
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The integration constants are h, h, a, h u Je lt a 2 . The symbols A and M 
stand for the following expressions: 

A =-2ffl^ 8 + 2(mc z + h)tf + 2{m l -\-m z )^ — 2(hc 2 — Jc)tf 

— 2(m 1 + m z )c z ?J i — (2k + a 2 + (3)c z X i + (3c i , 

M = — 2ni(i s + 2(mc 2 + h)ix a + 2(m 1 —m 2 )(i 5 —2(hc 2 — Jc)(i i 

— 2 (wij — m 2 ) c V s — (2 h + a 2 + /?) c> 2 + /? c* ; 

or, setting for short 

A = 2 m, 5 = 2 »? c 2 + h, C x — 2 {m x + w 2 ), C 2 = 2 (wjj — ?n 2 ), 

D = 2 (A c 2 — &), ^ = 2 (mj + «? 2 ) c 2 , ^ = 2 (^ — m 2 ) e 2 , 

^=2Jfe + a* + /8, G = (3c i 

A = — A%* + BV+ C 1 -k i —Dtf — E 1 tf — F??+ G, 

M = —A(i s + Bn e + C z (i 5 — Dfi i — F z[l 3 —Ffi 2 + G. 

It was shown by Serret* that under certain conditions 

AM =0 

is a singular solution of the problem of Euler or of that of Lagrange, and as his 
reasoning is directly applicable to the above problem we may conclude that this 
equation is also a singular solution of the problem in hand. In the original 
problem of Euler A and M are of the fourth degree in 1 and p respectively, in 
Lagrange's problem they are of the sixth, and in this problem they are of the 
eighth. 

The condition for such singular solution is that X — X , \ being the initial 
value of \ be a multiple factor of A of order m, m>2; or similarly, that ft — (i 0) 
/i being the initial value of ^, be a multiple factor of M of order n, n>2; that is, 
a, must be at least a double root of the equation A =0 ; or fi at least a double 
root of M = 0. 

The trajectory in the meridian plane will therefore be the ellipse a, = X , 
or the hyperbola (i = fx 0} according as 

'3 AN 
'o 
or 

«. = ». ©.=»• 

'8A^ 

'o 



A » = °> (§s). = °. 



where A and ( -*— ) are the values of A and its derivative with regard to % 

\OA/o 



*(Euvres de Lagrange, Vol. XII, Note II, pp. 346-349. 
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when a, = ;\, and likewise M and { -~— ) are the values of M and its derivative 

with regard to (x. when ^ = ^ . 

As the meridian plane rotates about the line of centres, the ellipse X = yl 
will generate an ellipsoid of revolution and the hyperbola (i = (i will generate 
an hyperboloid of revolution. Upon one of these surfaces the motion takes place. 

But if /L = /l and p = ^ at the same time where A and (i occur at least 
twice as roots of the equations A = and M = respectively, then no motion 
is possible in the meridian plane, for in this case 

d% d(i _ 

dt ~ ' dt ~ ' 

in space, however, the trajectory will be a circle whose axis is the axis of revo- 
lution and whose centre is K z . This kind of motion was first remarked by 
Legendre,* and later by Andrade.f 

With regard to equations of the form 

'dx\* 



(m) = *<•> 



Charlier J established certain theorems which are immediately available for our 
discussion. The equation is supposed to be one capable of interpretation in a 
mechanical problem, and therefore x as well as its first derivative with regard to t 
may be regarded as real, finite, and continuous for all finite values of t. 

In the first place let a be any real root whatever of F{x) = of order m, 
m being a positive integer. Let x = x at the beginning of motion and let x 
be so taken that no zero value of F{x) lies between a and a; . During the motion 
let x approach a continuously by increasing or decreasing values according as x 
is less than or greater than a. Under these assumptions, Charlier succeeded in 
establishing the following results : 

First, for w > 2, x can not reach the value a in a finite interval of time, 
and being real and continuous, it can not therefore pass the value a. 

Second, for m = 1, x will reach the value a within a finite time ; but it can 
not pass this value, for then the derivative will change sign, which is not 
admissible. 

* TraiU de fonctions elliptiques, Vol. II, p. 539. 
t Journal d'Ecole Poly technique, Cahier 60, 1890. 
% Meehanik des Simmels, Vol. I, pp. 117-163. 
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In other words, x can never pass a zero value of F(x). 

In the second place, let a and b be two consecutive roots of F (x) = of 
orders rn and n respectively, let a < b and let x lie between a and b at the 
beginning of motion. Then it follows immediately from the above theorem 
that if x = x at the beginning of motion, it must always remain within the 
interval (a, b). If m = n = 1 and if x starts from a- toward a, it will continue 
to decrease until it reaches the value a ; then it turns and continues to increase 
until it takes the value 6, when it changes its direction again toward a, and so on 
indefinitely. The motion in this case is one of oscillation, or what Charlier calls 
libration motion. The limits a and b are called the Kbration limits. 

If m>2 and n = 1, and if x starts from x , x being in the interval (a, b) 
at the beginning of motion, toward a, it will continue to decrease but will never 
reach the value a ; it approaches a as a limit and can not attain this limit in a 
finite interval of time. If x starts from x toward 6, it will reach J in a finite 
interval of time ; it will then change its direction and continue to decrease 
toward a as before. If m>2 and n>2, then, in whichever direction x moves 
at the beginning of motion, it will continue to move indefinitely, approaching 
either a or & as a limit. Charlier calls this kind of motion limitation motion and 
the limits, limitation limits. 

These theorems will now be applied to equations (30) and (31). 

The coordinates A, and p have been introduced by the equations 

where A, and (i are the parameters of a system of confocal conies having their 
real foci at the fixed centres of force, K lf K 2 . If we denote the distance K x K z 
by 2 c, then the least value of 1 is + c> while the greatest and least values of fx 
are +c and — c respectively. Since A, is the semi-major axis of the ellipse, 
the semi-minor axis is V/L 2 — c 3 ; of the hyperbola the semi-transverse axis is p, 
and the semi-conjugate axis is Vc 2 - (i z . Hence the equations of the ellipses 
and hyperbolas are respectively 



±i 4. y — i jm — y — i 

J^2 s i8 „2 *> ..2 „2 ..2 — ± > 

and therefore 



^ T A a — c 2 — ' a 2 c 2 — << 2 



H u o" — (i* 
c* x 2 = tf (i* , c 2 y 2 = (A 2 — c 2 ) (c 2 — p 2 ) . 
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From these defining equations of X and (i, it follows that 9i and fi must 
satisfy the inequalities 

00>/1>C, +C>(M> — c. 

In order that the quantities under the radicals in the integrals may be 

positive, that is, in order that there may be real motion, X and (i must also 

satisfy the conditions 

A > , M > , 

respectively. All of the above conditions % and p must satisfy throughout 
the motion. 

The potential function TJ can be written 

p =nra{(?(^)- 2<sv+4) )-(p(^)- 2 < s '' 2 +*))}- 

It is evident that, of the possible values /I and (i, TJ becomes infinite for 
"k = c and for ^ = ± c or zero. Therefore "k and (i can not be allowed to take 
these values nor to approach them indefinitely, if we wish to avoid infinite dis- 
continuities in TJ. Hence there must be a lower limit to the possible values 
of /I, as l^>c; likewise, an upper and a lower limit to the possible positive 
values of (i which we denote by \ and l{ respectively, and an upper and lower 
limit to the possible negative values of (i which will be denoted by V z and ^ 
respectively. The following inequalities must therefore be fulfilled by JL and p 
respectively : 

°o>a>7, h>(*>l'i, V % >ti>\. 

For any value of % >-c for which A vanishes, the function TJ will be finite, 
however little this value of % may differ from c. Hence the least value of the 
possible values of % will always be less than or at most equal to the least value 
of X for which A vanishes. In other words, all the positive roots of the equation 
A = which are greater than + c are also greater than or equal to I. Similarly, 
all the positive values of (i in the interval (0, +c) for which M vanishes lie in 
the interval (l[, ZJ, and all the negative values of (i in the interval ( — c, 0) for 
which M vanishes lie in the interval (? 2 , ?|). 

When n becomes zero, not only TJ but also d(i/dt becomes infinite. But for 
this value of (i, no discontinuities occur in the three final integrals of the solution ; 
that is, the time t as well as the trajectory in the meridian plane and in space 
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is finite and continuous, as one can easily determine. Hence in a discussion of 
the trajectories no further attention need be given to these discontinuities in TJ 
and d/i/dt and it will be sufficient that % and (i satisfy the conditions 

cc ~^> %>l f l 2 > /* > Zj . 

Without loss of generality we may take in 1 >■ m % . Then, since c 2 and a 2 
are always positive, if we regard m, m 1 , m 2 and /3 as essentially positive, no 
change of sign can take place in any of the coefficients of A and M which are 
independent of the arbitrary constants h and h. Hence A, C lf C % , E lt E % , Gr 
are all positive, while B, D, F may be positive or negative, depending for their 
sign upon the values of h and h. 

When /l = 0, the expression A is positive, and when 2= ± c, A is negative. 
A must therefore have one real root between and + c, and one between — c 
and 0. By the aid of Descartes' Rule of Signs, we find also that A can not have 
more than three positive, and it can not have more than five negative roots. 
A = must therefore have at least one positive root, but may have three; and it 
must have at least one negative root, but may have three or five. The same is 
true relative to the roots of the equation M = 0. 

Further, for a certain value of %, the function A becomes negative and 
remains negative for all greater values of A. Hence there must be a greatest 
value of % for which A is positive or zero ; that is, there must be a greatest value 
for which A changes from positive to negative. This value, which we shall 
denote by a x , must be a root of A = 0. Let us suppose a x > c. Since A is 
negative for %■=■ + c and positive for values of % less than a lf A must also 
change sign for some value of 2 between + c and a x ; there must be a positive 
root of A=0 in the interval (+c, a a ). It is evident then, that A is positive 
for all values of % in the interval (a 2 , %) and for no other values of 2 greater 
than + c. Since a x and a z are both greater than Z, X may take all values in the 
interval (a a , a x ) including the limits a 2 and a x . 

We have seen that one positive root of A = lies between and + c. If 
a x be this root, then there can be but one positive real root ; if a x be different 
from this root but less than + c, then a s will also be less than + c. In neither 
of these cases is motion possible. In all cases of real motion, therefore, we must 
have both a a and a 2 greater than + c, and X must satisfy the inequality 
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As to the roots of M = we have remarked that those only can have any 
influence upon the motion which satisfy the inequality l^^n^li, al) d that M = 
must have one positive root but may have three, and one negative root but may 
have three or five. We denote the positive roots by b u b Z) b z and the negative 
roots by (3 lt (3 2 , @ 3 , fl if /3 6 . By b t we shall always mean the single positive 
root, or, if there be three positive roots, the greatest of the positive roots in the 
interval (0, Z a ) for which M changes from positive to negative values. Similarly, 
by /?! we mean the negative root which must always exist, or, if there be more 
than one negative root, the greatest numerically of the negative roots in the 
interval (4, 0) which causes M to change from negative to positive values for 
increasing values of (i. 




Fio. 5. 

Hence it is clear that the trajectory of the moving body M in the meridian 
plane must remain wholly within one of two regions enclosed by Jl = a lf h = a 2 
and by the hyperbola branches (i-=(3 lt ^.z=.b 1} as shown in Fig. 5. As the 
meridian plane rotates about the line of centres, these regions will generate a 
tube-like portion of space within which the trajectory must always remain. 

Considering the definitions and restrictions with regard to the values of (i 
for which M vanishes as determined above, it is found that there are 105 different 
combinations of the roots of the equation M = which must be taken into 
account in a complete discussion of the motion, and each of these combinations 
must be associated with each of the two combinations of the roots of the equation 
A = 0, which are . ^ . __ 

thus making in all 210 different combinations of the roots in question, which 
may influence the motion. 
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The positive roots of M = may occur in any one of the following ways : 
i?, : b % , b 3 imaginary, or real and greater than l lt 



J i 
A 
A 
A 
A 



b z , b$ real, unequal and less than b x , 

b 2 , b s real, equal and less than b 1} 

b 2 , b 3 real, equal, greater than J x and less than \, 

b %> b s real and equal to b t . 



The ways in which the negative roots may occur relative to /^ and \ are 
accounted for in the following manner: 

First, any two of them, as /? 2 , @ 3 : 



G 6 



/3 2 , /? 3 imaginary, or real and less than 4, 

$5, /? 3 real, unequal and greater than (3 1} 

/? 2 , /? 3 real, equal and greater than (j lf 

/2 2 , @ 3 real, equal, less than /3 X and greater than \ , 

fi % , (3 3 real and equal to fix . 



Second, the remaining two, /? 4 , /? B : 



A 
A 
A 
A 
A 



/? 4 , /? B imaginary, or real and less than \ , 

/? 4 , /3 B real, unequal and greater than $ x , 

/3 4 , /? 6 real, equal and greater than ^ , 

/? 4 , /? 6 real, equal, less than /? a and greater than ? 2 , 

ft » ft rea, l an ^ equal to /?j . 



We must also take certain of the (7's with certain of the As in order 
to obtain relations between the negative roots which are not given above. 
These are : 

real, unequal, greater than /?] , and /? 2 = /? 3 > /? 5 > /? 4 , 
real, unequal, greater than /3 X , and /3 6 > /3 3 = /3 3 >■ /3 4 , 
real, unequal, greater than /?! , and /9 6 > /? 4 > /? 3 = /? 2 > /3 X , 
real, unequal, greater than ^ , and /? 6 = /? 3 = /? 2 > /? 4 , 
real, unequal, greater than & , and /? 6 > /? 4 = /? 3 = /? 2 , 
real, equal, greater than /3 X , and /? 2 = /? 3 = /? 4 = /? B , 
real, equal, greater than /Sj , and /? 6 = /? 4 > /3 3 = /? 2 > ^ , 
real, equal, less than (3 lt greater than \ and /3 B = /3 4 =/3 3 = /3 2 , 
real, equal, less than (3 lt greater than ? 2 and ft = ft>ft = ft. 



^3 Al 


ft, ft 


C 3 D ZZ 


ft, ft 


C^D^ 


ft, ft 


@z A* 


ft, ft 


^3 Ad 


ft, ft 


^3 Al 


ft, ft 


Cb Aa 


ft, ft 


CiAi 


ft, ft 


4 A. 


ft, ft 
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All the combinations of the roots of A = and M = 0, which are actually 
different and which can have any influence on the motion, can now be formed. 
They are given by the following formulae : 

Pa k i = -o-i -Bj C k D x , P ij3 (2j k) = A % Bj G 3 Z/ 2i k , 

Pijm2 = -^-i Bj C m D 2 , Pijp (j>, l) = A % Bj C p D Pt j , 

Pijn3 == A i B j C n D S) Pij&n == •A- t B J C 6 D nf 

where 

i = 1, 2 ; / = 1, 2, 3, 4, 5 ; h = 1, 2, 3, 4, 5 ; 

? = 1, 2 ; «i = 2, 4, 5 ; rc = 4, 5 ; p = 3, 4 . 

These different combinations, 210 in all, give rise to numerous possibilities of 
motion, depending upon the initial values of % and (i ; but, as one can easily 
determine, all the possible cases of motion so arising will fall into one of the 
following six general classes : 

I. Libration in Pi and libration in (i, 

II. Libration in % and limitation in ^, 

III. Libration in Pi and (i = constant, 

IV. Pi = constant and libration in fi, 
V. % = constant and limitation in [i, 

VI. Pi = constant and (i = constant. 

To give a clear idea of the nature of the motion in the various possible cases 
it will be sufficient to discuss a few of them so as to illustrate each of the above 
six general classes. To this end we shall consider further only the combinations 

-^133(2,1) aD d P 233(2,1)- 

In the first of these combinations we have 

«i>« 2 , b 3 = b 2 <b u & = &>&>&, 
and therefore the following possibilities of motion are presented : 

(1) a x > \ > a 2 : b x > <u > b 2 = b 3 ; 

(2) a 1 >?l >a 2 : fi =b 2 = b 3 ; 

(3) a x > X > a 2 : b 2 = b 3 > (i > (3 2 = 0, ; 

(4) «! > \ > a 2 : p = (3 2 = fi 3 ; 

(5) a 1 >^ >a 2 : & = & > ^ > /? 6 ; 

(6) <*! > X > a 2 : Pi>l*o>Pi, 

Pl and ^ being the initial values of Jl and /* respectively. In each of the six 
47 
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cases there is libration in k. In (6) there is also libration in (i and therefore (6) 
belongs to class I. In (1), (3), (5) there is limitation in p as in class II, and 
(2) and (4) belong to class III. 

If at the beginning of the motion the conditions (6) are satisfied, the trajec- 
tory of the moving point in the meridian plane must always remain within one 
of the two trapezoidal regions bounded by the ellipses k = a lf k = a 2 and by the 
hyperbola branches fi = /? 4 , (t = /? a . Since there is libration in both k and (i, 
these variables will take the values a lf a s and /3 4 , /Jj respectively an indefinite 
number of times and the trajectory will meet the bounding arcs an indefinite 
number of times. As the meridian plane rotates about its axis, this region of 
the meridian plane will generate a tube-like portion of space, bounded by the 
surfaces of the ellipsoids of revolution k = a x , k = a 2 and the hyperboloids of 
revolution /x = /3 4 , [i = ^ t . Within this tube the trajectory must always remain, 
meeting the limiting surfaces an indefinite number of times. The trajectory 
will be a coil-like curve. 

If conditions (1) are satisfied at the beginning of motion, the trajectory 
in the meridian plane must always remain within the region bounded by the 
ellipses k = a lt k = a 2 and the hyperbolas (i = b lf (i = b s . The variable k will 
take the values a lf a z an indefinite number of times while p may take the 
value b t once, if at all, and must then approach b % as a limit, not reaching it in a 
finite interval of time. Hence the trajectory will meet the bounding arcs of the 
ellipses an indefinite number of times and approach the bounding arc of the 
hyperbola (i = b % in an asymptotic manner. In space the trajectory will be a 
wave-like, spiral-like curve meeting the surfaces of the ellipsoids of revolution 
k = a x , k = a % an indefinite number of times while it continues to approach 
asymptotically the surface of the byperboloid of revolution ^ = 6 2 , which is 
intercepted between the ellipsoids Jl = a x , "k = a g . 

When conditions (3) or (5) hold at the beginning of motion, the trajectory 
will be of the same general character as that corresponding to conditions (1). 

When (2) and (4) are satisfied at the beginning of motion, we have libration 
in "k and fj. = constant throughout the motion. In the meridian plane the trajec- 
tory consists of the arc of the hyperbola intercepted between the ellipses %, = a 1} 
k =a a and the motion is one of oscillation upon this arc. In space the trajectory 
must remain upon the surface of the hyperboloid p = b z intercepted between the 
ellipsoids generated by "k =a lf X -=.a % . It is therefore a wave-like curve crossing 
this surface and meeting the bounding circles an indefinite number of times. 
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In the second combination selected for further discussion, -P 23 3 (2)1 ), we have 
«i = «», h = h<h, Pz = Ps>p6>Pi>Pi, 

and therefore we have the following possibilities depending upon the initial 
values /l and fx of /I and [i respectively : 



(1) 


Ao = <i\ = a.% : 


&i > /"o > &a = h 5 


(2) 


^<o = a i == tt a • 


^o = h = ^3 ; 


(3) 


^o ~~ a i — "a • 


6 a = h > i"o > Ps = Ps ; 


(4) 


Aq — C£i — — (Xa * 


fo = Pa = P3 j 


(5) 


«0 = a l ■— a 2 ' 


Pn = Ps>Ho> Pal 


(«) 


Aq = flfj — 0! 2 : 


p5>Po> Pi- 



Of these different cases (6) belongs to class IV ; (1), (3), (5) to class V 
and (2), (4) to class VI. 

If conditions (6) are satisfied at the beginning of motion, the trajectory will 
consist of the arc of the ellipse % = a lt between the hyperbolas [i = /? B , (i = /? 4 . 
In space it will be upon the surface of the ellipsoid of revolution bounded by the 
circles of intersection of this surface with the two hyperboloids of revolution 
H = (3 5f fi = (3 i . The trajectory will be a wave-like curve upon the ellipsoid 
crossing its surface and meeting the bounding circles an indefinite number of 
times. 

If at the beginning of motion conditions (1) are fulfilled, the trajectory in 
the meridian plane will be the arc of the ellipse h = a x intercepted between the 
hyperbolas (i = b l , [i = b 2} but the motion will not be one of oscillation as in the 
preceding case. If the variable (i move initially in the direction of 6 a , it will 
reach that value in a finite time, then change its direction toward b z , which it 
approaches as a limit, not reaching it in a finite interval of time. As the meridian 
plane rotates about its axis, the motion will be upon the surface of the ellipsoid 
of revolution a, = a x , between the hyperboloids of revolution (i=zb 1} p = \ • 
The trajectory will be a curve winding about the ellipsoid and approaching 
asymptotically the circle of intersection of this surface with the hyperboloid 
H = b 2 . In cases (3) and (5) the trajectory is of the same nature as that just 
described. 
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In cases (2) and (4) no motion is possible in the rotating plane, but in space 
the trajectory is a circle whose axis is the line K X K Z . In (2), for example, 
the trajectory is the circle of intersection of the ellipsoid of revolution X = a x 
and the hyperboloid n = b 2 . 

By giving particular values to the masses and the arbitrary constants, 
interesting particular cases of motion may be discovered. If we take a = 
and make the necessary changes throughout, the cases of motion in a fixed plane 
passing through the line of centres will be obtained. Then, if we take >l = c 
throughout the motion, the trajectory will consist of the straight line segment 
R X K % ; or, if we take fi = ± c throughout the motion, the trajectory will 
consist of the segments of the line K x K% beyond the fixed centres K % and K x 
respectively. 

Pkincbton, N. J. 



